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Preface 


In  this  thesis,  techniques  are  presented  which  will  aid 
in  the  study  of  scattering  from  a  loaded  four  arm  planar 
equiangular  spiral  antenna.  Extensive  use  was  made  of  the 
references  listed  in  the  bibliography  and  the  serious 
reader  is  encouraged  to  study  those  references.  By  under¬ 
standing  the  content  of  the  references,  and  therefore  the 
structure  on  which  this  thesis  is  built,  the  reader  should 
be  able  to  extend  the  techniques  described  here  to  other 
equiangular  spiral  antennas  and  possibly  to  other  frequency 
independent  antennas. 

I  would  like  to  acknowledge  my  thesis  advisor.  Dr.  A,  J. 
Terzuoli  for  extending  to  me  the  flexibility  required  in 
this  creative  effort.  I  wish  also  to  thank  my  wife,  Kelly, 
for  her  multifaceted  support. 


Stephen  Charles  Moraites 
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Abstract 


This  thesis  presents  analytical  and  numerical  techniques 
for  analyzing  the  scattering  from  a  planar  equiangular  spiral 
antenna  loaded  at  its  terminals.  The  scattering  matrix  for 
the  loaded  antenna  is  derived  as  a  function  of  the  antenna 
load  impedance.  This  derivation  is  the  result  of  an  analytical 
study  of  the  voltage/current  relationship  at  the  antenna 
terminals  along  with  an  application  of  a  multiport  analysis 
of  the  scattering  problem.  A  numerical  technique  is  also 
developed  that  utilizes  at  different  wavelengths  common 
elements  of  the  generalized  impedance  matrix  of  the  moment 
method  solution.  This  technique  provides  for  rapid  computation 
of  the  scattering  properties  of  the  spiral  antenna  across  a 
band  of  frequencies. 
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CHAPTER  1 


Introduction 


This  thesis  pertains  to  the  electromagnetic  scattering 
properties  of  a  loaded  four  arm  planar  equiangular  spiral 
antenna . 

The  spiral  antenna  is  one  of  a  class  of  antennas  with 
frequency  independent  properties  (17),  As  a  radiating 
device,  its  properties  have  been  measured  (A),  analyzed 
(17),  and  exploited  (9;21).  As  a  scatterer  however,  no 
published  work  has  yet  fully  characterized  this  structure. 
This  thesis  addresses  two  questions  on  the  problem  of 
determining  the  scattering  properties  of  the  loaded  four 
arm  planar  equiangular  spiral: 

1,  What  effect  does  the  antenna  load  have  on  the 
scattering  properties? 

2.  How  can  the  special  geometry  of  the  symmetric, 
self  complementary,  equiangular  spiral  be  used 
in  the  numerical  analysis  of  this  scatterer? 

These  two  questions  will  be  dealt  with,  respectively,  in  thi 

thesis  by  accomplishing  the  following: 

1.  An  expression  is  derived  which  gives  the 
scattering  matrix  of  the  antenna  as  a  function 
of  antenna  load  impedance.  The  expression  is 
in  terms  of  theoretically  determined  parameters 
and  numerically  determined  parameters  (from  the 
method  of  moments.) 

2.  Techniques  are  developed  that  decrease  the 
computation  time  of  the  numerically  determined 
parameters.  These  techniques  are  developed 
from  the  special  geometric  properties  of  the 
antenna . 


The  importance  of  antenna  scattering  theory  lies  in  its 
applicability  to  remote  prediction  of  antenna  parameters, 
field  measurements,  radar  cross  section  studies,  and  electro¬ 
magnetic  compatibility.  Previous  work  pertaining  to  antenna 
scattering  is  readily  found  in  the  literature.  Various 
approaches  to  studying  antenna  scattering  are 

1.  spherical  mode  expansion  of  the  incident  field 
(3:317-329) 

2.  multiport  systems  analysis  (6:107-125) 

3.  equivalent  circuit  analysis  (5:1-23) 

4.  receiver  antenna  current  analysis  (5:1-23) 

5.  numerical  computation;  e.g.  moment  method 
solution  of  the  field  equations  (6) 

Green's  work  (5)  provides  a  survey  of  approaches  1  through  5 

and  is  a  standard  reference  on  the  subject. 

This  thesis  approaches  the  problem  of  antenna  scattering 

by  describing  the  scattering  situation  as  a  multiport  system 

whose  parameters  are  found  both  analytically  and  through 

numerical  computation,  A  complete  analysis  of  the  multiport 

requires  a  knowledge  of  the  port  parameters  that  relate  the 

port  currents  to  the  port  voltages.  Chapter  2  of  this  thesis 

presents  the  derivation  of  the  port  parameters  for  the  four 

arm  planar  equiangular  spiral  antenna.  Also  within  Chapter  2 

is  information  on  the  radiation  properties  of  the  antenna 

that  is  pertinent  to  the  subsequent  scattering  analysis. 

Chapter  3  formulates  the  problem  of  scattering  by  a  loaded 

four  arm  spiral.  Here  the  measurement  port  is  added  to  the 

multiport  system  and  the  re.Ta '  cnship  among  the  various  port 


parameters  is  explained.  The  end  result  of  Chapter  3  is  an 
expression  for  the  scattering  matrix  as  a  function  of  the 
scattering  antenna  load  impedance.  This  expression  is  in 
terms  of  scattering  matrices  of  the  antenna  with  its  terminals 
shorted  in  various  ways.  Each  of  these  "shorted"  scattering 
matrices  can  be  computed  numerically  using  the  method  of 
moments . 

Chapter  4  briefly  explains  the  moment  method  as  applied 
to  the  spiral  antenna.  Techniques  are  developed  that  exploit 
the  very  special  geometry  of  the  planar  equiangular  spiral 
and  allow  one  to  compute  the  scattering  properties  of  the 
antenna  over  a  band  of  frequencies  with  a  minimum  of  compu¬ 
tation  time. 

Chapter  5  is  the  conclusion  and  includes  suggestions  on 
the  implementation  of  the  analytical  tools  provided  by  this 
thesis. 
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CHAPTER  2 


The  Spiral  Antenna 


Geometry 

A  four  arm  planar  equiangular  spiral  antenna  is  drawn 
in  Figure  2.1.  The  spirals  forming  the  edges  of  the  arms 
are  described  by  the  polar  equation  2.1  (A:182;  17:14) 

0=0o  +  [  tan(A)  ]llnp  (2.1) 

or 

p=a  exp[a(0-6)]  (2.2) 

where  p  and  0  are  standard  polar  coordinates  and  the  remain¬ 
ing  constants  are  parameters  that  determine  the  expansion 
rate  and  angular  position  of  the  spiral.  From  the  form  of 
equation  (2.1),  it  is  apparent  that  two  angles,  0o  and  A 
completely  define  the  geometry  of  an  infinite  spiral.  The 
performance  of  an  antenna  whose  geometry  is  entirely  deter¬ 
mined  by  angles  is  independent  of  frequency  (17:13).  If  the 
arms  of  an  infinite  spiral  antenna  are  truncated  so  as  to  be 
physically  realizable,  the  spiral  structure  still  performs 
--  in  an  antenna  sense  —  as  the  infinite  structure  does, 
provided  that  the  frequency  of  operation  remains  within 


certain  bounds. 


Bandwidth 


Frequency  independent  performance  of  the  infinite  spiral 
antenna  means  1)  the  antenna  terminal  impedances  are  indepen¬ 
dent  of  frequency;  and  2)  the  radiation  pattern  is  independent 
of  frequency  (17:13).  Over  the  bandwidth  of  the  truncated 
spiral  antenna,  the  antenna's  impedance  and  pattern  will  be 
treated  as  independent  of  frequency.  Using  the  frequency 
bounds  stated  by  Mosko  (9:98),  the  upper  limit  is  given  by 
X^=TTdi  and  the  lower  limit  by  Xj^=Tid2/3  where  dj  and  d^are  as 
shown  in  Figure  2.1.  (These  bounds  permit  multimode  operation) 
Thus  the  antenna  shown  in  Figure  2.1  has  a  6  to  1  frequency 
range . 

Terminal  Impedances 

Before  considering  terminal  impedances  of  a  four  terminal 
spiral  antenna,  one  should  be  familiar  with  the  impedance  rela¬ 
tionship  of  single  port  planar  complementary  antennas.  (Planar 
complementary  antennas  are  a  pair  of  structures  that  when 
superimposed,  one  on  the  other,  exactly  fill  a  plane  with  no 
overlap.)  Given  two  such  antennas: 

ZZ^=n|/4  =  (60TT)2  (2.3) 

where  Z  and  Z^  are  the  input  impedances  of  an  antenna  and  its 
complement  (2;S371),  and  no  is  the  characteristic  impedance 
of  free  space.  If  an  antenna  is  its  own  complement,  then 
Z  =  Z^  sbOiT.  Deschamps  extended  this  relationship  to  multi¬ 
terminal,  self  complementary  structures  (2).  Referring  to 


Figure  2.2,  which  shows  the  voltage  and  current  conven¬ 
tions  for  a  four  terminal  spiral  antenna,  the  relationship 
between  voltages  and  currents  is,  from  Deschamps, 


or 
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[u]=[c]Cl] 


(2.5 


Application  of  Deschamps*  work  to  a  four  arm  planar  equi¬ 


angular  spiral  antenna  requires  evaluation  of  the  following 
transformation: 


Co 
Cl 
C2 
C3 

where 

=ri  0 /4sin(mTT /4)  m=l,2,3 

5o  =0 

The  reader  should  refer  to  reference  2  for  details  of  the 
derivation  of  equations  (2.4)  and  (2,6),  The  evaluation  of 
equation  2.6  yields 
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From  a  mathmatical  point  of  view,  the  system  of  equations  2,4 
has  meaning  only  on  a  hyperplane  imbedded  in  four-space.  In 
the  derivation  of  equation  2.4,  the  requirement  was  made  that 
and  ^(,^=0.  Therefore,  there  are  only  three  independent 
coordinates  of  the  vector  space  spanned  by  [u]or  [l]  .  It 
will  prove  useful  to  apply  a  transformation  of  coordinates 
that  allows  one  to  arbitrarily  fix  the  value  of  one  transformed 
component  of  [  i.  ]  to  some  convenient  value.  The  transformation 
is  a  simple  rotation  in  four-space  about  the  origin.  Its 
physical  interpretation  is  shown  in  Figure  2.3.  The  four 
terminal  currents  are  represented  by  the  vector  [1].  Now 
define  branch  currents  [I]  such  that 


To"]  r  1  0  0  -ri  rt" 

li  -1  1  0  0  Ii 
i.2  “  0-110  I2 

I3J  L  0  0  “1  Ll  L^i 


(2.8) 
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Terminals  numbered  clockwise  O  fhru  J 


? — I 


—Vo¬ 


'S  ^ 


-V^-4^ 


Figure  2.3.  Branch  Current  and  Branch  Voltage  Conventions. 


(2.9) 


And  since  the  matrix  relating  [  i,  ]  to  [I]  is  of  rank  three, 
the  complete  solution  for  [I]  can  be  expressed  as  (9:211,219) 


1  0  +  1 1 
to  +  n  +  1 2 
0 


(2.10) 


which  shows  that  [lalmay  be  arbitrarily  fixed. 

To  simplify  the  subsequent  analysis,  the  voltage  vector 
[uj  will  be  transformed  also.  A  new  vector  as  shown  in 


Figure  2.3  is  defined  by 


v  1  - 

V2 

V, 


1-10  0 
0  1-10 
0  0  1-1 
10  0  1 


(2.11) 


[V]  =  [Trtul 

Equations  (2.5),  (2.9)  and  (2.12)  can  be  rewritten  as 


where 


Evaluating  [Z] 


[V]=[Zl[I] 


[Z]=[Tftc][T] 


(2.12) 


(2.13) 


[Z]  = 


Zq  Zj  Z2  Z3 
Z3  Zfl  Z I  Z  2 
Z  2  Z3  Z(J  Zy 
Z I  Z  2  Z3  Zo_ 


(2.15 


where 


2Co-Cx-C3  1  +  * 

2Ci-C2-Co  „  -1 

2  C2-C3-C1  ~  1-. 

2C3-C0-C2  -1 


(2.16: 


The  physical  system  described  by  equation  2.13  is  shown  in 
Figure  2.3. 


Radiation  Modes 

By  proper  phasing  of  the  terminal  currents,  [  I  ]  ,  one 
can  obtain  radiation  patterns  that  are  suitable  for  direction 
finding  systems  (9).  Selection  of  three  orthogonal  phase 
schemes  for  the  terminal  currents  leads  to  three  independent 
modes  of  radiation.  Two  of  the  modes  are  shown  in  Figure  2.4 
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MODE  1 


• 

1 

-J 

Figure  2.4.  Radiation  Patterns  of  Four  Arm  Spiral  in  Plane 
Perpendicular  to  Plane  of  Antenna.  Phase  of 
Source*  (9-92)^®^"’^’'^^^  Currents  Shown  to  Right  of  Patterns. 

Feed  Point  Configuration 

To  obtain  the  two  modes  shown  in  Figure  2.4,  one  must 
choose  a  method  of  feeding  the  antenna.  Either  mode  could 
be  created  by  feeding  the  structure  with  two  balanced  trans¬ 
mission  lines.  Figure  2.5  shows  a  feed  method  for  either 
mode.  This  method,  however,  cannot  be  used  if  both  modes 
are  simultaneously  required.  If  both  modes  are  required, 
than  the  feed  method  shown  in  Figure  2.6  can  be  used. 
Physically,  this  method  can  be  realized  by  running  a  coax 
cable  to  each  terminal  and  connecting  the  outer  conductors 
together.  The  impedance,  Z,  in  the  figures  represents  the 
internal  impedance  of  the  source  as  seen  from  the  antenna. 
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Figure  2.5 


Antenna  Feed  Networks  tor  Single  Mode 
Operation . 


When  receiving  energy  from  space,  the  antenna  also  sees  an 
impedance,  Z,  caused  by  the  receiver  input  impedance  and 
the  transmission  line  between  the  receiver  and  the  antenna. 
Throughout  this  thesis,  the  impedance  seen  at  the  antenna  by 
looking  back  through  the  transmission  line  and  into  any 
receiver  circuitry  is  called  the  antenna  load  impedance. 
Chapter  3  is  devoted  to  the  effect  of  this  load  impedance  on 
the  antenna's  scattering  properties.  The  input  impedance  for 
a  multiterminal  antenna  depends  on  the  mode  being  radiated 
(9:100).  From  reference  9,  or  by  use  of  equation  2. A,  mode  1 
antenna  impedance  is  133  ohms.  Mode  2  antenna  impedance  is 
9A  ohms.  This  is  the  impedance  seen  by  each  generator  (and 
generator  impedance)  in  Figure  2.6, 

Radiation  Pattern 

In  free  space,  a  planar  spiral  antenna  radiates  on  both 
sides  of  the  plane  of  the  antenna.  This  thesis,  however, 
deals  with  an  antenna  mounted  on  the  surface  of  an  absorbing 
cavity.  Such  antennas  are  constructed  to  eliminate  reflec¬ 
tions  from  structures  behind  the  antenna.  The  polarization 
is  elliptic  such  that  the  electric  vector  rotates  as  the 
spiral  expands  (17:A0). 


CHAPTER  3 


The  Loaded  Spiral  Scatterer 

The  theory  of  scattering  from  a  loaded  four  arm  equi¬ 
angular  spiral  antenna  will  be  developed  from  more  general 
cases  of  electromagnetic  scattering.  The  formulation  of  the 
theory  follows  Harrington  (8)  but  is  modified  so  as  to 
pertain  to  backseat tering  only. 

Scattering,  No  Load 

Consider  a  measurement  antenna  and  the  environment 
around  it  as  a  single  port  network  shown  in  Figure  3.1.  By 
removing  any  scatterer  of  interest  from  the  environment  and 

relating  V_  to  a  specific  ratio,  Z _  is  obtained.  Under 

scatterless  conditions: 


V^/I 


r 


=Zo 

rr 


(3.1) 


Where  the  superscript  indicates  no  scatterer  of  interest 


is  the  input  impedance  of  the  measurement  antenna.  Now 
if  a  scatterer  is  placed  in  the  environment,  the  voltage 
and  the  current  are  related  by 


V 


s 

r 


/I 


(3. 


Within  the  network  of  Figure  3.1,  a  wave  is  radiated  from 
the  measurement  antenna,  impinges  on  the  scatterer  and  is 
scattered,  with  some  of  the  scattered  energy  appearing  as 
a  voltage, AV^  ,  at  the  terminals  of  the  measurement  antenna: 

AVj.  =V®-V“  (3. 


A  measure  of  the  scatterer's  effectiveness  in  reflecting 
energy  back  to  its  source  is  AZ^,  the  scattering  impedance. 

AZr  =AVj./Ij.  (3,. 


The  familiar  measure,  radar  cross  section, g 
AZ^  as  follows  (7:167) 


a  =Att 


2 


is  related  to 


(3. 


where 

R  =  distance  from  scatterer  to  measurement  antenna 
rio=  characteristic  impedance  of  free  space 
i,  =  effective  length  of  measurement  antenna 
X  =  wavelength  of  incident  energy 


In  the  analysis  to  follow,  o  is  an  inadequate  measure  of  the 
scatterer  since  phase  information  is  lost  in  its  definition. 
Note  that  the  development  so  far  also  assumes  some  fixed 


polarization  for  both  the  incident  and  scattered  energy, 
the  two  being  not  necessarily  the  same.  A  more  complete 
description  of  the  scatterer  is  a  scattering  matrix  (8:619) 
which  relates  the  incident  field  to  the  scattered  field  E®; 


JYR  exp(jkr) 
i)fR  exp(jkr) 


S 11  S 12  El 

S21  S22  E2 


where  the  subscripts  refer  to  two  orthogonal  polarizations. 

S  and  AZ^  are  related  by 
mn  r  •' 

=,nn=  ^  (3. 

Note  that  the  dependence  of  AZ^  on  m  and  n;  that  is,  on  the 

polarization  of  interest,  is  not  explicit  but  is  just  as 
definite  as  in  the  S  parameters. 


Loaded  Single  Port  Scattering 

Now,  returning  to  the  network  representation  of  the 
scattering  problem,  one  can  conceptualize  the  antenna  scat¬ 
tering  problem  by  adding  an  additional  port  to  represent  the 
terminals  of  the  scatterer  as  in  Figure  3.2. 


Figure  3.2.  Single  Port  Scatterer  Added  to  Network  of  Fig3,l 


A  matrix  equation  relating  the  port  parameters  is  (7:166) 


(3.8) 


where  AV^  is  defined  in  equation  3.3  and  the  Z  variables 

are  standard  open  circuit  impedance  parameters.  If  the 

scatterer  is  loaded  with  z.  at  its  port,  then  v  =-Z  I  snd 

L  s  L  s 

equation  3.8  can  be  solved  for  the  scattering  impedance, 


AZ  (Z.  )  : 
r  L 


AZ^(Z^^)=  ^ 


AZ. 


rr 


z 

Z 

rs 

sr 

Z 

+  Zt 

ss 

L 

(3.9) 


Loaded  Two  Port  Scattering 

If  the  scatterer  has  two  ports,  the  relationship  among 
the  port  parameters  is  (8:618) 
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the  load  relations, 
impedance : 


(3.10) 


and  solving 


where 

Z^^=Z^  =  the  open  circuit  impedance  parameter  between 

the  measurement  antenna  and  port  i 
^ij  ”  circuit  impedance  parameter  between  port  i  and  j 

7^  .  =  load  impedance  at  port  i 


Figure  3.3.  Network  Representation  of  Scattering  by  a  Two 
Port  Scatterer  Loaded  with  and 

Equations  3.10  and  3.11  can  be  used  to  represent  the  loaded 
four  arm  spiral  antenna  only  if  a  definite  terminal-to-port 
assignment  is  made  and  only  if  no  self  complementary  two  port 
theory  (2:S377)  is  applied.  (Because  all  the  terminals  meet 
at  a  point,  the  feed  configuration  at  one  port,  however  it's 
defined,  will  affect  the  feed  configuration  of  the  other  port 
when  the  complement  of  the  antenna  is  considered.  The  four 
terminal  equiangular  spiral  is  not  a  self  complementary  two 
pert  structure.) 

Multiterminal  Loaded  Scattering 


As  explained  in  Chapter  2,  a  feed  system  more  versatile 
than  the  two  port  may  be  required.  The  network  representing 
the  four  terminal  scattering  system  is  shown  in  Figure  3.4. 


Figure  3.4.  Network  representation  of  Scattering  by  a  Four 

Terminal  Loaded  Structure.  (Measurement  port  on  left) 


Making  use  of  equation  2.13,  one  can  arrive  at  the  system  of 
equations  describing  Figure  3.4. 


(3.12) 


or  by  separating  the  system  and  using  the  load  constraints: 


Parameter  Study 

Before  proceeding  further  to  express  equations  3.13  and 
3.14  into  more  practical  forms,  one  can  gain  insight  into 
antenna  scattering  by  considering  the  meaning  of  the  various 
terms  in  the  preceeding  equations.  With  reference  to  equation 
3.8,  the  parameter  Z  is  related  to  the  scattering  antenna 
power  gain  G  (8:620): 


(3.15) 


where  the  constant  of  proportionality  includes  measurement 
antenna  parameters,  distance,  etc.  so  that  G  is  a  character¬ 
istic  only  of  the  scattering  antenna's  radiation  properties. 
Another  parameter,  AZ^^.  is  the  value  of  AZ^.  when 


(3.16) 


In  the  analysis  to  follow,  Z^^  and  its  multiport  counterparts 
are  unknown  quantities  and  it  is  convenient  to  omit  them  from 
the  equation.  (Using  the  method  of  moments,  described  in 
Chapter  4,  calculation  of  Z  requires  the  introduction  of 
generators  and  wires  in  the  antenna  model.  This  addition  is 
undesirable.)  Solving  for  Z^^  in  terms  of  open  and  short 
circuited  scattering  impedance  yields 


Z2 

rs 


=Zgs(AZr(<»)  -  AZ^(O)) 


(3.17) 


Substitution  and  further  manipulation  lead  to  expressions 
similar  to  those  derived  by  Green  (5:24-26): 


iZ_.(ZL)=- 


Z,iZ,(»)  ♦  Z„AZ,(0) 

Li  L  SSL 


z  +  z, 

ss  L 


AZ  (Z, )=AZ  (Z  *)  + 
r  L  r  ss 


Zs,[AZ^(»)  -  AZ^(0)](Zl  -  Z^*) 


(3.19 

(3.2C 


(Z  +  Z  *)(Z  +  Z.  ) 

ss  ss  ss  L 


Equation  3.20  allows  one  to  define  a  structural  and  antenna 
component  of  the  scattered  field.  AZ  (Z  *)  of  the  right 
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hand  side  is  termed  the  structural  component  and  the  remainder 
of  the  right  hand  side  is  termed  the  antenna  component  since 
it  goes  to  zero  when  the  antenna  load  impedance  Zj^  equals  the 
conjugate  of  the  antenna  input  impedance,  Z^^. 

Referring  now  to  equation  3.10  describing  the  two  port 
representation,  Z^^  and  Z^^  are  related  to  the  radiation 
properties  from  their  respective  ports  of  the  scattering 
antenna.  When  the  voltage  at  the  measurement  antenna  is 


so 


V  =I  Z 
r  in 


(3.21 


G=C| I.Z  +  I  Z  1  2  ^3  22 

I  j  ri  2  ra  '  vj.zz 

Again,  solving  for  Z^^  in  terms  of  open  and  short  circuited 
scattering  impedances: 


Zi.  =  Z,,  rAZ„(oo,oo)  -  AZ_(a>,on 


(3.23 


and 


(3.24) 


therefore 


G=C 


VZj  1  UZj.(-oo)  -  AZj.(0,<»)lIj  -  AZj.(«,0)]I^j  2 

Further  manipulation  of  the  two  port  equations  yields 
’^Lz ^  ~ 

^Li^La  Z(«,oo)+ZgZL2AZ(0,«>)+ZgZL^AZ(»,0)  +  (Z|-Z2)AZ(0,0)  ^,^26) 

where 


(' 


Zg=Zjj=Z22  of  equation  3.10 
equation  3.16 

Note  the  assumption  that  the  ports  are  defined  such  that 

shown  from  equation  2.13 

that  the  two  port  impedances  are  real.  Therefore,  a  conjugate 

load  impedance  is  real  and  equal  to  some  Z^.  Thus 

Zj  Z(oo,0)  +  AZ(0,co)]  +  ZjAZ(oo,«)+  Z(0,0)] 

AZ  (Z.,Z.)=-2 - 1 -  (3.27) 

r  i’  1  2(Z  +Z  ) 

^  s  i' 


which  is  the  structural  component  of  the  scattered  field  for 
the  particular  terminal  to  port  assignment.  Note  that,  unlike 
a  single  port  scatterer,  a  two  port  scatterer  does  not  have 
a  unique  structural  scattering  component  since  the  terminal 
configuration  and  so  the  conjugate  load,  Z^,  can  vary. 


22 


.25) 


Referring  finally  to  equations  3.13  and  3.14  describing 
the  more  general  loaded  four  terminal  structure,  the  transfer 
impedances  are  analogous  to  their  two  port  counterparts: 


Z2  =Z  [AZ  « 
r  0  0  r  oo 


-  AZ. 


(3.28) 


ri  0  rl«> 


A  Z^  8  ) 
r  oo 


(3.29) 


Z"  =ZJAZ  « 

r  2  0  r  oo 


r,  OO  , 

AZ  0  ] 

r  0 


(3.30) 


Z2  =z  [AZ  « 
ra  0  r  “ 


(3.31) 


G=C  I  Z,  +  I  Z,  +  I  Z  +  I  Z  ^ 
0  ro  in  2  r2  a  ra 


(3.32) 


In  order  to  express  equation  3.13  in  terms  not  containing 
Z^^,  one  can  make  use  of  the  results  of  Chapter  2  and  let 
1 3=0.  Also,  the  last  equation  of  the  system  of  equations 
3.14  can  be  dropped  since  V3=-(Vo +Vi +V2 ) .  Appendix  A 
contains  the  details  of  the  algebra.  The  following  results 
apply  to  the  case  where  Zj^^=Zj^  i=0,l,2,3  and  use  is  made  of 
the  equality  between  Zj  and  Z3  of  equation  3.12: 


where 

Fi  =  Zo(Zl+(Zo+Zi ))2/D=Zo/(Zl+Zi ) 

Fa  =  2(Zo+Zi)[(ZL-Zi)2-Zn/D 
D  =  [Zj^+(Zo+Zi)]2[Zl-Zi] 

Equation  3.33  expresses  the  scattering  impedance,  AZ^,  as  a 
function  of  the  antenna  load  impedance,  Z^.  The  expression 
is  in  terms  of  analytically  determined  parameters  Z^  from 
equation  2.16  and  numerically  determined  parameters  A Z^ [] . 
Since,  as  will  be  explained  in  Chapter  4,  the  scattering 
matrix  of  equation  3.6  is  the  representation  of  the  scatterer 
most  easily  obtained,  equation  3.33  can  be  rewritten  in  scat¬ 
tering  matrix  notation.  From  equation  3.7,  each  element  of 
S  is  related  to  its  corresponding  A  Z^  by  the  same  factor. 


As  a  result,  the  matrix  S  is  a  function  of  the  antenna  load 
impedance  and  is  equal  to 


where  Fi  and  F2  are  as  in  equation  3.33. 

Regarding  the  concept  of  structural  and  antenna  components 
of  the  scattered  field,  each  radiation  mode  will  have  its  own 
input  impedance  (9:100)  therefore,  the  proper  conjugate 
matched  load  will  depend  on  the  mode  under  consideration. 

Since  a  single  impedance  value  for  cannot  provide  a  match 
for  all  modes  simultaneously,  the  concept  of  a  structural 
component  does  not  have  the  meaning  it  had  for  the  single 


CHAPTER  4 


Field  Computation 
Formulation  of  the  Problem 

The  end  objective  of  the  field  computation  at  hand  is  to 
have  a  numerical  value  for  a  complex  matrix  S  (equation  3.6) 
that  expresses  the  magnitude  and  phase  of  the  scattered 
electric  field  at  a  point  in  space.  The  problem  of  obtaining 
this  end  objective  is  divided  into  three  subproblems: 

1.  Select  the  incident  field 

a.  wavelength 

b.  polarization 

2.  Find  the  induced  surface  current  on  the  scatterer 

3.  Compute  the  electric  field  due  to  the  surface 
current 

Then,  by  making  use  of  equation  3.6,  the  scattering  matrix, 

S,  can  be  computed.  The  numerical  technique  used  to  accomplish 
the  necessary  computation  is  the  moment  method.  Throughout 
this  thesis,  Newman's  Electromagnetic  Surface  Patch  (ESP) 
code  (12)  is  the  implementation  of  the  moment  method  that  is 
used  as  a  reference.  If  one  makes  use  of  this  code  directly, 
then  subproblem  1  is  normally  taken  care  of  by  the  required 
inputs  to  the  code.  The  solutions  to  subproblems  2  and  3 
and  the  computation  of  S  are  standard  outputs.  If  one  needs 
to  make  use  of  the  techniques  explained  in  this  chapter, 
however,  then  some  reorganization  of  this  code  is  required. 

For  this  reason,  the  three  subproblems  are  identified  and 
explained  seperately.  Since  the  rational  for  choosing  a 


method  of  solving  each  subproblem  is  based  on  the  subsequent 
step,  the  subproblems  will  be  explained  in  reverse  order. 


Computation  of  the  Fields 

From  field  equivalence  theorems  (18:518,519),  the 
electromagnetic  field  reflected  from  the  surface,  S,  of  a 
material  body  is  equivalent  to  the  field  generated  by  appro¬ 
priate  surface  currents  on  S  in  free  space.  In  the  case  of 
a  thin,  flat,  perfect  conductor,  the  equivalent  surface 
currents  on  the  two  flat  sides  may  be  considered  equal  and 
coincident.  A  single  electric  surface  current  on  the  open 
surface  of  a  zero  thickness  plane  can  be  used  to  find  the 
reflected  fields  if  the  single  current  is  the  vector  sum 
of  the  two  surface  currents  of  the  flat  plate  (11:784),  The 

electric  field  from  a  surface  current  J  (7)  is  given  by 

s 


E(7) 


r(r,r ’  )J 


s 


(7*  )ds* 


L(J^(r')) 

S 


(4.1) 


where 

r_  is  the  free  space  Green's  function 

r  is  symbolic  for  the  complete  set  of  spatial  coordinates 
L  is  a  linear  operator  relating  surface  current  to  the 
electric  field 

Jg  is  assumed  known  as  an  expansion  of  subsectional  basis 
functions  (6:11): 


where  each  is  defined  only  over  a  small  patch  of  the  entire 
surface  and  is  the  complex  amplitude.  If  is  chosen  as 
a  sinusoidal  spatial  variation  in  the  current  density,  then 
equation  (A.l)  is  known  in  closed  form  ( 1 5 ; 18 : 370 , 401 ) . 
References  10  and  11  explain  the  use  of  patches  of  sinusoidal 
surface  currents  to  model  surfaces.  The  approach  to  finding 
the  total  scattered  electric  field  by  solving  for  fields  of 
sinusoidal  surface  currents  on  small  patches  of  the  scatterer 
requires  that  the  surface  to  be  modeled  be  subdivided  into 
patches.  Appendix  B  explains  how  to  subdivide  the  spiral 
antenna. 

Given,  then, the  functional  form  of  and  its  amplitude, 
the  fields  from  all  surface  patch  currents  are  calculated 
using  the  appropriate  formula  as  in  equation  (4,1),  The 
total  scattered  field  is  just  the  vector  sum  of  the  fields 
generated  by  1  =  1  to  M,  where  M  is  the  number  of  surface 

patch  currents. 

Computation  of  Surface  Currents 

Reaction.  The  induced  surface  currents  on  a  scatterer 
are  found  by  solving  an  integral  equation.  In  this  thesis, 
the  Electric  Field  Integral  Equation  (EFIE),  is  considered. 

The  EFIE  is  a  specific  form  of  the  Reaction  Integral  Equation 
(RIE)  (12:6).  The  moment  method  computer  code  used  as  a 
reference  (12)  makes  use  of  the  RIE. 

Reaction  is  a  measure  of  the  effect  of  one  electromagnetic 


source  upon  another  (16).  In  the  special  but  pertinent  case 


of  only  electric  surface  currents  existing  on  a  scatterer, 
the  reaction  between  two  sources  "a"  and  "b"  is 


-//^a  • 


(4.3) 


where  J  is  the  source  (surface)  current  on  "a"  and  E,  is 
a  b 

the  field  at  "a"  from  a  source  at  ”b“.  Reaction  is  reciprocal 
(16:1483): 


<J,,E.  >  =  <Jk,E  > 
a  b  b  a 


(4.4) 


Moment  Method.  From  Maxwell's  equations,  it  is  known  that 
the  electric  field  tangential  to  a  perfect  conductor  is  zero. 
Specifically,  the  tangential  electric  field  on  the  surface 
of  a  spiral  antenna  is  zero.  If  the  electric  field  around  an 
antenna  is  decomposed  into  an  incident  and  scattered  field; 
that  is,  if 


E  =  E^  +  E® 


then  the  boundary  conditions  imply 


E,.  =  E^  +  E®  =  0 

tan  tan  tan 


(4,5) 


(4,6) 
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Consider  E  to  be  radiated  from  some  unknown  surface 
current  on  the  antenna.  Using  the  linear  operator  of 
equation  (A.l): 


L(Jg)  =  E® 


(A. 8) 


Now  applying  equation  A. 8  to  equation  A. 7 


=  -E" 


;an 


tan 


(A. 9) 


The  method  of  moments  converts  equation  A. 9  into  a  system 
of  equations  whose  solution  is  an  approximation  to  (6). 


a; 

Let 

measure 

JU 

J  and 

» 

s 

between 

s  * 


As  a 


aid  of  equation  (A. 9), 


/A 


■  -A. 


•  E^  ds 
tan 


(A. 10) 


As  explained  in  the  beginning  of  this  chapter  by  equation 

A. 2,  Jg  is  expressed  as  a  sum  of  subsectional  basis  functions; 


(A. 11) 


i=l 


Substituting  this  into  equation  A. 10, 
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(A. 12) 
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where  the  "tangential"  notation  of  equation  4.10  is  dropped 
since  the  dot  product  of  the  electric  field  with  the  test 
current  (which  is  defined  to  exist  only  on  the  surface  of 
the  antenna)  selects  only  the  tangential  component.  Now  M 
linearly  independent  test  sources,  i=l  to  M  can  be  defined 
so  that  M  linearly  independent  equations  of  the  form  of 
equation  A. 12  result.  With  the  aid  of  the  reaction  notation, 
this  system  is 

"<Wj,L(Jj)>  <Wj,L(J2)>  •  •  •  <W^ ,L(Jm)>1  ri  1  r<-W^,Ei>" 

<W2,L(Jj)>  <W^,L(J2)>  <-W^,Ei> 


_<W^,L(J^)> 


<W„,L(J„)>J  LImJ  L<-WM»Ei>. 


[Zg][l]  =  [E] 

The  matrix  in  equation  A. 13  is  called  the  generalized  impe¬ 
dance  matrix;  its  elements  are  generalized  impedances  (6:8A). 
The  remainder  of  this  chapter  contains  special  techniques  for 
constructing  and  solving  equation  A. 13  specifically  for  the 
planar  equiangular  spiral  structure. 


Antenna  Symmetry.  As  shown  in  Appendix  B,  the  number  of 
modes  needed  to  represent  the  open  circuited  spiral  structure 
is  AN  where  N  is  the  number  of  modes  per  arm.  The  generalized 
impedance  matrix  contains  an  element  for  every  combination  of 
expansion  and  test  modes.  Therefore,  since  the  number  of  test 


modes  equals  the  number  of  expansion  modes,  there  are  (AN)^ 
elements  in  the  matrix.  By  making  use  of  the  antenna  symmetry 
however,  one  can  reduce  the  number  of  unique  elements  to  A(N)^ 
While  the  generalized  impedance  matrix  can  be  arranged  arbi¬ 
trarily,  i.e.,  the  modes  can  be  numbered  arbitrarily,  consider 
the  matrix  shown  in  Figure  A.l.  There  the  matrix  is  ordered 
and  partitioned  to  emphasize  the  arm  to  arm  impedances. 
Submatrix  contains  the  N*  generalized  impedances  between 

the  test  modes  on  the  i^.^^  arm  and  the  expansion  modes  on  the 
arm.  Because  of  the  rotational  symmetry  of  the  antenna, 
the  submatrices  on  the  same  diagonal  are  equal.  Thus  there 
are  A(N)^  unique  elements. 

Matrix  Symmetry.  If  the  test  modes  are  identical  to  the 
expansion  modes  in  both  functional  form  and  location,  then 

<W^,L(Jj)>  =  <J^,L(W^)>  (A.IA) 

and  since  reaction  is  reciprocal, 

<X,L(Wj)>  =  <Wj,L(X)>  (A. 15) 

from  which  it  follows  that 

<W^,L(Jj)>  =  <X,L(X)>  (A. 16) 

For  a  matrix  with  no  special  properties,  symmetry  obviously 
implies  that  the  number  of  unique  elements  in  the  generalized 
impedance  matrix  is  half  that  of  a  matrix  without  symmetry. 

For  the  matrix  represented  in  Figure  A.l  however,  the 
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Note: Matrix  symmetry  is  contingent  on  test /expansion  mode  equality. 


Figure  4.1.  Generalized  Impedance  Matrix  bhowing  Effect  of 
Antenna  Symmetry  (by  partitioning)  and  Matrix 
Symmetry  (by  shading). 

effect  of  symmetry  is  not  obvious  since  the  rotational  symmetry 
of  the  antenna  has  eliminated  the  need  to  compute  any  element 
below  the  diagonal  —  except  for  submatrix  An  .  Though  not 
obvious,  matrix  symmetry  does  further  reduce  the  number  of 
unique  elements  by  a  factor  of  two.  The  explanation  is  aided 
by  Figure  4.1,  where  one  can  confirm  that  the  elements  of  the 
submatraces  [A],  [B],  [C],  and  [D]  obey  the  following  equalities 


Sid 


-s 


a  .  .=a  . . 
b.  .=d  .  . 
c  .  .  =c  . . 
d  .  .=b  .  . 


Thus  if  the  upper  right  triangle  of  each  submatrix  is  computed, 
the  lower  left  triangles  are  obtained  from  the  above  set  of 
equalities.  Although  the  number  of  unique  elements  is  decreased 
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by  choosing  the  test  modes  to  be  the  same  as  the  expansion 
modes,  the  computation  time  for  each  unique  element  may 
increase  and  outweigh  any  advantages  offered  by  the  symmetric 
matrix.  As  explained  in  Appendix  B,  the  moment  method  computer 
code  used  as  a  reference  (12)  represents  the  modes  by  filaments 
of  current.  The  reaction,  the  generalized  impedance,  between 
two  modes  is  the  weighted  sum  of  the  reactions  among  the 
individual  filaments.  If  the  expansion  modes  are  represented 
by  K  filaments  and  the  test  modes  are  identical  to  the  expan¬ 
sion  modes,  then  there  are  filament  to  filament  reactions 
that  require  computation  for  each  unique  matrix  element.  If, 
however,  the  test  modes  are  represented  by  a  single  filament, 
then  there  are  only  K  reactions  to  compute  for  each  unique 
matrix  element.  So,  while  choosing  the  test  modes  to  be 
identical  with  the  expansion  modes  will  decrease  the  matrix 
computation  time  by  a  factor  of  two  by  virtue  of  reciprocity 
and  symmetry,  it  will  increase  the  time  by  a  factor  of  K  by 
virtue  of  the  number  of  test  mode  filaments.  The  net 
increase,  then,  of  time  is  a  factor  of  K/2  over  the  single 
filament  test  mode  case. 

Terminal  Short  Circuits.  Chapter  3  explained  the  need  for 
the  computation  of  the  scattered  field  with  the  feed  network 
configured  as  a  short  circuit  between  various  combinations  of 
two  terminals.  The  model  for  a  short  circuit  between  terminals 
is  explained  in  Appendix  B,  The  modeling  technique  requires 
the  addition  of  one  mode  for  each  short.  The  generalized 


impedance  matrix  will  be  altered  by  the  addition  of  a  row 
and  a  column  containing  the  reactions  between  the  shorting 
mode  and  the  modes  of  the  open  circuited  antenna  plus  one 
element  for  the  mode's  self  reaction.  The  generalized 
impedance  matrix  for  the  short  circuited  antenna  is  shown 
in  Figure  4,2,  Assuming  that  the  matrix  equation, 

[Z  j[l]=[Ei],  for  the  open  circuited  antenna  was  solved,  the 

O 

new  equation  can  be  solved  with  relatively  few  additional 
calculations.  Let 


[lol  =  [Z„l"  [El] 


(4.17) 


represent  the  solution  to  the  open  circuited  antenna  scat¬ 
tering  problem,  where  [Z  ]~^[Ei]  is  symbolic  of  the  solution, 

S 

regardless  of  the  method  used  to  obtain  it.  If  a  new  row 
and  column  are  added  along  with  an  additional  unknown,  I2, 
and  an  element,  Ea,  in  the  right  hand  side  vector,  then  the 
new  system  is  as  shown  in  Figure  4.2,  The  corresponding 


equations  are 


[Zg][lil  +  [Blla  =  [Ell 

[C][Ii]  +  DIa  =  Ea 


(4.18a) 


(4.18b) 


Solving  for  Ig  and  [Ii]  yields 


I,  = 


Ea  -  [C][li] 


(4.19) 


[Ill  =  [Z„]'‘[Ei]  -  [Z„]‘‘[B]I: 


(4.20) 
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Figure  A, 2,  Generalized  Impedance  Matrix  of  Spiral  Antenna  with  one 
Shorting  Mode. 


Applying  equation  A, 17  to  equation  A, 20,  substituting  into 
equation  A. 19  and  solving  again  for  I2  yields 


E2  -  [Cl[lJ 

I  - - 

D  -  [C][Z 

o 


and  substituting  into  equation  A. 20: 


-1 

[I  ]  =  [I  ]  -  [Z  1  ‘[B]  - 2 - - 

§  D  -  [C][Zg]  [B] 


(A. 21 


(A. 22 


If  the  solution  to  the  "short  circuit"  equation  is  obtained 
without  use  of  equations  A. 21  and  A. 22,  the  total  number  of 
numerical  computer  operations  is,  assuming  LU  decomposition 
and  M  modes,  approximately  1/3  M^  (1:152-160).  Making  use 
of  equations  A. 21  and  A. 22,  one  can  obtain  the  "short  circuit 


solution  with  approximately  2M*  operations.  Note  that  the 
evaluation  of  [Z  ]  * [Bl  can  be  made  along  with  other  open 
circuit  calculations. 

Selection  of  the  Incident  Field 

Wavelength.  The  scattering  properties  of  the  planar 
equiangular  spiral  structure  can  best  be  understood  if  its 
behavior  over  a  range  of  wavelengths  is  known.  By  making 
use  of  the  special  geometry  of  the  planar  equiangular 
spiral  one  can  avoid  the  need  to  recompute  the  entire 
generalized  impedance  matrix  of  equation  4.13  at  each 
different  wavelength.  The  structure  shown  in  Figure  2.1 
represents  an  antenna  with  a  6  to  1  wavelength  of  operation 
range.  Modeling  of  this  structure  through  the  use  of 
equation  3.33  must  be  done  within  this  range  since  the 
assumptions  about  frequency  independent  input  impedances  are 
not  valid  outside  this  range.  Also,  one  must  insure  that  the 
sizes  of  the  monopole  patches  are  no  greater  than  one  fourth 
wavelength  (See  Appendix  B).  Within  this  band  of  wavelengths, 
then,  the  following  explanation  is  valid. 

For  the  purpose  of  explanation,  the  modes  of  the  open 
circuit  spiral  structure  are  grouped  into  sections,  each 
section  comprised  of  the  dipole  modes  on  each  arm  correspond¬ 
ing  to  the  same  radial  distance.  For  instance,  the  innermost 
four  modes  —  one  on  each  arm  —  are  in  section  1.  The  next 
four  modes  farther  out  along  the  arms  form  section  2,  and  so 
on.  (This  assumes  that  only  the  current  modes  in  the  direction 


Figure  A. 3.  Fateh  Representation  of  Spiral  Showing 
Assignment  of  Patches  into  Sections  for 
Purpose  of  Illustrating  the  Scaling  of 
the  r.orieralized  Impedance  Elements. 


of  the  spiral  are  of  interest.  If  a  transverse  component 
is  also  modeled,  then  each  section  will  contain  an  additional 
four  modes.)  Figure  A. 3  illustrates  the  assignment  of  modes 
into  sections.  Although  two  surface  patches  are  associated 
with  each  dipole,  the  section  number  is  placed  in  the  patch 
in  which  the  dipole  begins.  With  this  sectioning  of  modes 
the  generalized  impedance  matrix  can  be  partitioned: 


where  the  superscripts  of  Z  refer  to  the  arm  number  as  shown 
in  Figure  A. 3.  The  argument  X  emphasizes  that  the  generalized 
impedances  are  computed  for  a  specific  wavelength  of  incident 
radiation.  The  impedance  entries  in  equation  A.2A  are  sums 
of  dipole  filament  to  dipole  filament  impedances,  each  of 
which  is  a  sum  of  monopole  filament  to  monopole  filament 
impedances  (20:13).  A  closed  form  expression  exists  for 
each  such  impedance  (lA).  The  exact  expression  is  not  of 
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interest  here  but  its  form  is.  In  a  lossless  medium,  the 
generalized  impedance  between  two  monopoles  can  be  expressed 
as  a  function  of  two  types  of  arguments: 

Z  =  Z(|,4')  (4.25) 

both  of  which  describe  the  geometry  of  the  two  monopoles. 

T  is  symbolic  for  several  variables  all  of  which  are  linear 
measures,  and  y  is  the  angle  between  the  two  monopoles 
(See  Figure  C.l.)  Now  for  any  x, 

Z(|,y)  =  Z(-^,4')  (4.26) 

That  is,  the  generalized  impedance  of  two  monopoles  at  some 
Xi  is  equal  to  the  generalized  impedance  of  two  monopoles 
scaled  by  x  at  a  wavelength  xXi,  assuming  f  remains  constant. 
Appendix  C  contains  the  details  of  the  proof  that  the  geometric 
relationship  between  two  monopoles  is  simply  a  scaled  version 
of  the  relationship  between  the  adjacent  set  of  monopoles. 

The  scaling  factor,  K,  for  adjacent  sets  of  monopoles  is  a 
function  of  the  angular  increment.  A,  between  adjacent  monopoles 

K  =  exp(aA)  (4.27) 

where  a  is  defined  in  equation  2.2.  Thus,  noting  that  the 
linear  dimensions  of  the  monopoles  are  implicit  in  the 
subscripts  of  equation  4.24,  the  elements  of  equation  4.24  and 
the  submatrices  of  equation  4.23  obey  the  following  equalities: 


'(i  +  l)(j  +  l) 


VN.A  1 ; 


and 


A .  .  ( X 1  )  - 


(A. 29) 


The  relationship  between  two  generalized  impedance  matrices 
for  two  wavelengths  related  by  a  factor  K  is  shown  below. 

If  the  generalized  impedance  matrix  is 


(A. 30) 


then 


(A. 31) 


where  [A]  contains  the  generalized  impedances  between  (N-1) 
sections.  Thus,  given  that  the  matrix  of  equation  A. 30  is 
known,  the  matrix  of  equation  A. 31  requires  only  the  addi¬ 
tional  computation  of  [P],  [Q]  and  [R] . 


Polarization.  It  is  known  from  experimental  work  (A) 
that  the  polarization  of  a  wave  radiated  from  a  spiral  antenna 


is  elliptical  with  an  axial  ratio  ranging  from  0  to  1.  This 
fact  suggest  that  the  scattering  analysis  should  proceed 
using  two  orthogonal  elliptically  polarized  waves  as  the 
components  of  the  incident  and  relected  waves.  Reinforcing 
this  idea  is  work  by  Wang  (21)  which  showed  the  utility  of 
expressing  the  scattering  properties  of  a  spiral  structure 
in  terms  of  left  and  right  hand  circularly  polarized  waves. 
In  Wang's  work,  the  phase  of  the  scattered  field  was  easily 
predicted  only  if  the  fields  were  expressed  as  circularly 
polarized  waves.  The  moment  method  computer  code  that  is 
used  as  a  reference  for  this  thesis  (12),  however,  expresses 
the  field  in  terms  of  0  and  $  components  of  a  spherical 
coordinate  system.  If  we  view  the  scattering  matrix,  ISj^], 
in  the  0,  4>  coordinate  system  as  a  linear  operator,  then  the 
scattering  matrix,  [S^]  ,  in  the  (G+j'I'),  (O-jl*)  coordinate 
system  is  easily  found  to  be 

[S^]  =  [B]"‘[Sl][B]  (4.32) 

where  the  columns  of  [B]  are  the  vectors  (l,j)  and  (1,-j), 
that  is,  the  new  basis  vectors  (9:332).  Therefore  the 
scattering  matrix  for  a  circularly  polarized  wave  is 


where 

S,  .  .  =  scattering  matrix  element  of  equation  3.6  for  a 
ui  J 

linearly  polarized  (0,$)  wave. 


CHAPTER  5 


Conclusions 

Techniques  for  analyzing  the  scattering  properties  of  a 
loaded  four  arm  planar  equiangular  spiral  antenna  were  pre¬ 
sented.  The  two  main  results  of  this  thesis  are  the  deriva¬ 
tion  of  the  scattering  equation  for  the  loaded  antenna  in 
Chapter  3  and  the  development  of  a  technique  to  make  use  of 
the  moment  method  generalized  impedance  elements  at  different 
wavelengths  in  Chapter  4.  Implementing  these  techniques 
requires  a  working  knowledge  of  a  moment  method  code  like 
that  of  reference  12  and  some  careful  "bookkeeping". 

Topic  for  Further  Work 

The  result  of  Chapter  4,  a  technique  that  makes  use  of 
the  same  generalized  impedances  at  different  wavelengths,  was 
developed  for  a  self  complementary  structure.  A  possible 
extension  of  this  work  to  large  planar  sections  could  be 
developed  by  dividing  a  plane  section  into  an  equiangular 
spiral  and  its  complement  and  then  applying  the  technique. 
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Appendix  A 

Derivation  of  the  Scattering  Impedance 

The  objective  of  this  appendix  is  to  derive  an  expression 
for  AZ^  in  terms  of  load  impedances,  port  parameters,  and 
short  circuited  scattering  impedances.  The  expression  will 
be  developed  from  equation  3.12,  repeated  here 
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(3.12) 


Solving  for  AV^.  yields 


AV  =  AZ  I  +Z  I  +Z  I  +Z  I  +Z  I 
r  rr  r  ro  o  ri  i  rz  z  rs  j 


(A.l) 


and  for  AZ^,  the  scattering  impedance 

AV  f  Z  .1  . 


AZ. 


(A. 2) 


If,  in  equation  A.l,  I^  is  set  to  zero  and  the  other  currents 


I.  are  set  equal  to  one  another;  1^=1,  then  AV  =I^Z  . .  Fr 
equation  2.8,  implies  that  all  terminal  currents  are 

zero.  Therefore,  AV^,  which  is  a  measure  of  the  signal 
received  from  the  spiral,  is  zero  and 

f  Z  .=0 
k  ri 


om 


(A. 3) 
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Referring  now  to  equation  A. 2,  aZ^  is  a  function  of  and 
therefore  of  the  load  impedances.  This  functional  dependence 
will  be  denoted  by 


(A. A) 


Let  for  all  i.  Then  all  terminal  currents  are  zero  and 
from  equation  2.10,  the  branch  currents  for  all  i.  From 
equation  A. 2  and  A. 3  then, 


Now  let  Ij=I^  =  I^=0  in  equation  A.?,.  From  equation  2,8, 
terminals  2  and  3  must  be  open  circuited  .  If, 

in  addition  to  this  constraint,  Z.  *0  then  V,=0  and  from 

the  second  line  of  equation  3.12,  1 ^--Z ^^1 ^IZ Substituting 


In  a  similar  manner,  it  can  be  shown  that 


The  load  constraints  on  the  branch  voltages  are 
(See  Figure  3.4) 
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V3 

) 

+ 

Zl. <1.-13) 

(A. 13) 

Making  use  of  the  load  conditions,  letting  13=0  (see 
equation  2.10),  and  eliminating  the  last  equation  of  the 
system  in  equation  3.12  (since  M ^=-(V g+V 2) )  allows  one  to 
•S  express  the  system  as 
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Solving  equation  A. 14  for  [I]  and  substituting  into  equation 
A, 2  yields 


ro 


[Z  A  +Z  A  +Z  A  ] 
or  11  irzi  2r  3 1 ' 
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(A. 15) 
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A.  .  is  the  cofactor  of  the  iith  element  of  the  matrix 
i  i 

in  equation  A. 14 

A  is  the  determinent  of  the  matrix  equation  A.l4 


Z  .=Z. 
ri  ir 


One  can  now  solve  for  Z^^  in  terms  of  Z^()  with  open  and 
short  circuited  loads  at  the  terminals. 


(A. 16) 


(A. 17) 


and  since  Z  = -(Z  +Z  +Z  ) ,  equation  A. 9  becomes 
3r  or  ir  zr 


Substituting  equations  A. 6,  A. 7,  A. 8,  A. 16,  A. 17,  and  A. 18 
into  equation  A. 15  results  in  the  scattering  impedance  as 
expressed  in  equation  3.33. 


Appendix  B 

Surface  Patch  Modeling 

Functional  Form 

The  moment  method  code  used  as  a  reference  (12)  represents 
surface  current  as  a  sum  of  piecewise  sinusoidal  surface 
monopoles.  Two  adjacent  monopoles  form  a  dipole  and,  as 
illustrated  in  Figure  B.l,  a  reference  direction  is  assigned 
to  the  current  flow.  The  moment  method  provides  a  solution 
for  the  phase  and  magnitude  of  each  dipole.  Figure  B.l  shows 
six  rectangular  sinusoidal  surface  monopoles  joined  together 
to  form  three  overlapping  surface  dipoles.  By  adjusting  the 
magnitude  and  phase  of  each  dipole,  one  can  model  a  continuous 
current.  Each  dipole  is  called  an  expansion  mode  since  the 
current  is  "expanded"  into  a  linear  combination  of  these 
modes.  Non-rectangular  surface  currents  are  modeled  similarly, 
as  shown  in  Figure  B.2  (10).  Note  that  in  both  Figures  B.l 
and  B.2,  the  current  goes  to  zero  at  the  ends  of  the  dipoles 
so  that  no  accumulation  of  charge  needs  to  be  accounted  for. 
Though  not  shown  in  the  figures,  there  is  generally  also  a 
surface  current  running  transverse  to  that  shown  in  the 
figures. 

Numerical  Treatment 

If  the  surface  monopoles  are  exact  rectangles  then  the 
electromagnetic  fields  associated  with  them  have  closed  form, 
exact  expressions  (15).  Otherwise,  the  fields  are  found  by 


integrating  the  fields  from  filamentary  sinusoidal  monopoles, 
the  fields  from  five  filamentary  monopoles  are  summed  (12:27). 
As  expected,  the  number  of  integration  intervals  will  affect 
the  field  computation  time.  Work  by  Tulyathan  (20:13)  showed 
that  an  integration  interval  of  0.07  times  a  wavelength  is 
adequate  for  representing  a  surface  current.  Additional 
precision  is  obtained  only  with  an  increase  in  computation 
time.  For  the  antenna  structure  shown  in  Figure  2..  1,  four 
intervals  (five  filaments)  per  segment  is  sufficient  for 
modeling  the  component  of  current  in  the  direction  of  the 
spiral.  The  length  of  the  monopoles  also  affects  both  the 
precision  and  time  of  the  field  calculation.  Electromagnetic 
modeling  with  monopoles  as  large  as  one  fourth  wavelength 
yields  good  results  (20:37).  The  largest  patch  then  on  the 
spiral  on  which  a  monopole  is  superimposed  should  not  exceed 
one  fourth  wavelength.  For  the  spiral  antenna,  the  patch 
length  determines  the  number  of  monopole  patches.  The  antenna 
of  Figure  2.1  requires  N  patches  per  arm  or  4N  patches  total. 

Patch  Construction 

Two  methods  of  subdividing  the  spiral  structure  into 
monopole  patches  are  discussed  below. 

Skinner  constructed  the  arms  of  the  antenna  from  quadra- 
lateral  (trapezoidal)  patches  whose  corners  are  formed  by  the 
intersection  of  radial  lines  and  the  spiral  curves  (19:38). 
Figure  B.3  illustrates  this  method.  Note  that  the  angle 
between  adjacent  radials  determines  the  size  of  each  patch. 
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When  the  self  complementary  property  of  the  antenna  is  to  be 
modeled,  the  angle,  A  ,  must  be  a  submultiple  of  45  degrees 
so  that  a  rotation  of  45  degrees  brings  the  structure  back 
on  its  complement.  Although  division  of  the  arms  could  be 
made  at  different  angular  increments  (e.g.  45“  for  the 
inner  part  and  22^“  for  the  outer  part)  frequency  scaling 
of  the  generalized  impedance  elements  as  explained  in 
Chapter  4  is  more  straightforward  if  a  constant  A  is  used. 
If  A  =22^“,  then  the  antenna  of  Figure  2.1  can  be  modeled 
using  204  patches. 

Instead  of  using  straight  radials  to  subdivide  the 
spiral,  one  can  intersect  the  spiral  with  a  family  of 
orthogonal  curves.  The  equation  of  the  orthogonal  curve 
is  given  by  equation  B.l  , 

p=a“^a^  exp[-  ^(0  -  6)]  (B.l) 

where  all  parameters  and  variables  have  the  same  value 
and  meaning  as  in  equation  2.2.  Note  that  the  orthogonal 
curve  is  also  an  equiangular  spiral  but  it  expands  as  0 
becomes  more  negative. 

Test  Modes 

As  explained  in  Chapter  4,  the  testing  functions  can 
also  be  thought  of  as  surface  currents  in  a  reaction  with 
the  expansion  mode  currents.  For  these  testing  mode  currents 
one  can  choose  the  same  or  a  different  number  of  filamentary 
monopoles  to  represent  the  testing  modes  as  is  used  to 


represent  the  expansion  modes.  Regardless  of  the  number, 
the  functional  form  of  the  testing  currents  in  the  reference 
moment  method  code  (12:21)  is  a  piecewise  continuous  sinusoid 


Terminal  Shorting  Modes 

So  far,  explanation  of  the  expansion  and  testing  modes 
has  been  limited  to  those  modes  associated  with  the  open 
circuited  model  of  the  antenna.  An  additional  mode  is 
created  by  placing  a  monopole  on  one  arm  at  the  terminal  end 
of  the  spiral  and  another  monopole  on  the  other  arm  as 
illustrated  in  Figure  B.A.  The  two  monopoles  together  form 
a  single  dipole  mode  whose  magnitude  and  phase  can  be  solved 
for  as  explained  in  Chapter  4.  This  technique  simply  provides 
for  the  continuity  of  current  from  one  arm  to  another. 
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Figure  B.A,  Terminal  Ends  of  the  Spiral  Showing  the 

Construction  of  a  Shorting  Mode.  Distance 
between  Monopoles  must  be  Small, 


Figure  C.2  shows  two  monopoles,  I  and  J,  and  the  adjacent 
set,  I+l  and  J+1,  along  with  the  lengths  and  angle  corre¬ 
sponding  to  those  in  Figure  C.l.  To  support  the  text  in 
Chapter  4,  it  is  sufficient  to  show  that  triangles  AjOiBa 
and  B1O2A2  are  similar  triangles  and  that  the  sides  of 
A1O1B2  are  K  times  the  sides  of  B1O2A2  where  K  is  from 
equation  4.27.  The  similar  triangles  A2PB2,  B2PC2,  CjPBi 
and  BiPAj  imply  that  0i=02»  Therefore  Similar 

triangles  A2PB1  and  B2PA1  imply  that  angle  PA2B1  equals 
angle  PB2Ai  and  that  angle  AjB2C2  equals  angle  B1A2B2. 
Therefore,  triangles  A1O1B2  and  Bi02A2  are  similar.  And 
since,  by  the  equiangular  construction  (equation  2.2), 
PB2=KPA2,  then  the  equality  AiB2=KBjA2  is  obtained.  Thus 
the  scale  factor  of  the  two  similar  triangles  A1O1B2  and 
B1O2A2  is  K. 


Figure  C.2.  Adjacent  Sets  of  Monopoles  and  Associated 
Geometry. 
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